Draft version September 26, 2011 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



CRUSTAL FAILURE DURING BINARY INSPIRAL 



o 

(N 
D 

m 
(N 

00 
I 

o 



o 

O 



A.J. Penner^, N. Andersson^, D.I. Jones\ L Samuelsson^ & I. Hawke^ 

^School of Mathematics, University of Southampton, Southampton S017 IBJ, UK and 
■^Department of Physics, Umea University, SE-901 87 Umea, Sweden 
Draft version September 26, 2011 

ABSTRACT 

We present the first fully relativistic calculations of the crustal strain induced in a neutron star by a 
binary companion at the late stages of inspiral, employing realistic equations of state for the fluid core 
and the solid crust. We show that while the deep crust is likely to fail only shortly before coalescence, 
there is a large variation in elastic strain, with the outermost layers failing relatively early on in the 
inspiral. We discuss the significance of the results for both electromagnetic and gravitational-wave 
astronomy. 

Subject headings: Stars: neutron — Gravitational waves — dense matter 



1. INTRODUCTION 

Finite size effects for neutron stars in coalescing bi- 
nary systems are of great interest since the departure 
from point-mass dynamics may provide a unique insight 
into the high-density equation of state when such sys- 
tems a re detected by future gravit ational- wave experi- 
ments (jFlanagan fc Hindereii 120081) . In this paper we 
advance our understanding of this problem by providing 
the first general relativistic calculations of the strain in- 
duced in the solid crust during binary inspiral, building 
upon our recent rel ativistic elastic perturbation scheme 
(jPenner et al.l[20TTI ). 

The solid crust makes up a relatively small fraction 
of the star, and its shear modulus is small relative to 
the pressure. Nevertheless, there are several reasons why 
an accurate treatment of the crustal strain problem is 
of interest. Firstly, a sufficiently high signal-to-noise 
gravitational-wave observation may allow the effects of 
the solid phase to be identified. Secondly, large scale fail- 
ures of the crust during the late stages of inspiral may 
have interesting effects, especially given the link between 
such mergers and short gamma-ray bursts. Thirdly, the 
techniques developed here will prove essential in tackling 
the related problem of calculating, in General Relativity, 
the maximum mass quadrupole a compact star can sus- 
tain, a problem of great intere s t for gravitational- wave 
searches (see e.g. lAbbott et al.l (|2010[ )). 

Of particular interest is to know at what stage in the 
inspiral the crust fails. Early work by iKochaneM ()1992[ ) 
indicated this oc curs only shortly befo re r uerger. Re- 
cent estimates bv lPostnikov et all ()2010f ) and lTroia et al.l 
(|2010,) suggest failure at a much earlier stage. We re- 
solve this discrepancy below. However, our modelling 
indicates very different levels of strain in different parts 
of the crust, suggesting a much richer failure sequence 
than simple considerations would have suggested. This 
is a subtlety that requires accurate calculations of the 
kind described here. 

2. BACK-OF-THE-ENVELOPE ESTIMATES 

Before discussing relativistic strained stars, we will 
make a few simple estimates using Newtonian theory. 
As well as providing useful insights, these estimates will 
enable us to normalise the relativistic calculations that 



follow. 

Consider a star of mass M, radius R, a distance D from 
a binary companion of mass M; 



comp- 



We focus on the 

situation where the star is axisymmetric, being deformed 
away from sphericity such that its quadrupole moment 
is changed by a fractional amount e. This is non-zero by 
virtue of the tidal field of the companion, and also (we 
assume) by the solid phase having a zero-strain shape 
eo 7^ 0. This is obviously highly idealised, as in general 
the change in shape induced by the crustal strain will not 
coincide with the line joining the two stars. However, by 
considering this simple case we can make some useful 
estimates, and make contact with other results in the 

literature. 

The star's energy can be written ()Jone£ll2002[ ) 



E{e) = Eo + E, 



self— gravity 



tidal 



E. 



clastic 



(1) 



where Eq denotes the energy of the equivalent fluid star 
in the absence of a companion, -Escif-gravity the pertur- 
bation in the star's self-gravity, E'tidai the energy per- 
turbation due to the tidal interaction, and inelastic is the 
elastic energy. Since the self-gravity is minimised for 
a spherical configuration, we have itself -gravity Ae^, 
where A ^ GM^/R. Assuming that the elastic energy 
is minimised when the star has a shape e = eo, we have 
£^Giastic ^ B{e — eo)^, where B is of order the Coulomb 
binding energy of the solid phase (B ^ J fl dV , where fl is 
the shear modulus) . The gravitational field of the binary 
partner is, of course, -GAfcomp/^'comp, where rcomp = 
is the centre of mass of the companion. Expanding this 
about the centre of mass of the star of mass M, defined 
by r = 0, leads to the tidal field 



*tidai(r) = -r' ^3 P [ 



3n,n 



(2) 



where Ui is a unit vector along the line separating the 
two stars. The corresponding energy perturbation is 



^^tidal ~ $tidai(eAf) 



(3) 



Inserting these expressions into equation ([l} , minimis- 
ing with respect to e and making use of the expectation 
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that B <^ A leads to 
B 

'^A 



1 -^'^comp / R 



2 M 



D 



(4) 



The first term, proportional to en, is what is nor- 
mally termed the 'mountain' (jAbbott et al.l[20Toh . and 
is sourced by the star's non-spherical reference shape, 
caused by whatever geological processes may have been 
acting on it. The second term, proportional to Mcomp, 
is sourced by the tidal field created by the companion 
star. As expected in a perturbative treatment, these two 
contributions to the total shape e have decoupled neatly; 
we will denote them Cmountain and etidai- 

The corresponding perturbations in the star's 
quadrupole moment are found by simply multiplying by 
the moment of inertia /. Separating the pieces we have: 



^mountain — ^^^^ 

gtidai--2^^l;^l I 



(5) 
(6) 



Given that it is meaningful to separate the shape 
changes/quadrupole moments, it must be possible to 
separate the strains induced in the solid phases. In 
the mountain case, the actual deformation away from 
sphericity, emountain, is much less than the zero strain 
shape eo, by a factor B /A ^ 10~^, so the strain induced 
is; 

^mountain ^ ^0 ^mountain ~ 6o ~ (A/i?)eniountain (7) 

In contrast, for the tidal deformation, the only strains 
induced are by the shape change etidai itself, so that; 



Utidal ~ Etidai ^ Qtidal/-^ 



(8) 



We see that, for a given perturbation in the quadrupole 
moment, a factor of A/B ~ 10'°^ more strain is induced 
in the solid phases in the mountain case compared to 
the tidal deformation case. Physically, one way of un- 
derstanding this is as follows. In the mountain case, the 
deformation is sourced by the strains themselves, and 
the weakness of electrostatic forces versus pressure means 
one needs a large strain (~ eg) to build a given moun- 
tain (~ emountain)- In the tidal case, the deformation 
is sourced by the tidal field, with the consequent shape 
change (~ etidai) being the only source of strain; the 
solidity of the crust plays almost no role — the solid is 
deformed in much the same way as if it were fluid. 

It is now straightforward to estimate at what point 
during inspiral the crust fails. Combining ([6]) and ([8]) 
with Kepler's third law, i.e., il^ — GMtoted/ , where 
Q = 27r /orbit = T^fcw, we find that the crust fails 

(Utidal = "break) whcu 

106 cm /?^broak\l/2 



M 



1AM, 







1/2 



R 



0.1 



(9) 

(we have set Mcomp = M for simplicity). We have 
scaled the result, which should only be trusted up to 
factors of order unity, to the estimated breaking strain 
"break ~ 0.1 f rom rcccnt molecular dynamics simula- 
tions ([Horowitz fc Kadau. ,2009. ). Comparing with the 



estimated gr avitational- wave fre quency at the innermost 
stable orbit ()Abadie et al.ll2010D : 



-•isco 

JGW 



^63/2GMtotal 



= 1.6kHz(ii^) (10) 



(again, for an equal-mass binary) wc sec that crust failure 
should not be expected significantly before coalescence, 

in agreement with the results of Kochanek (1992). 

Thi s con trasts with the est imates of IPostnikov et al.l 
(|20T0l) and iTroia et all (poTol ) who found failure would 
occur much earlier in the inspiral. By using the map- 
ping between quadrupole deformation and strain appro- 
priate for mountains (effectively equation [T]) rather than 
that for tidal deformations (as per equation |5]), they 
overestimate the strain by a factor ^ A/B, altering 
the gravitational-wave frequency at failure by a factor 

~ (A/B)l/2 ^ 102. 

3. THE TIDAL PROBLEM 

The above estimates provide a useful guide, but they 
do not allow us to answer important questions such as 
where the crust first fails, and how large the strain varia- 
tions are within the crust. To a ddress these question s, we 
have extended the formalism of iPenner et al.l (1201 1[ ). As 
in that study, we take as our background model a spher- 
ical unstrained star (corresponding to eo = above) and 
subject it to a quadrupolar field, generated by a binary 
companion. However, in the present context we need to 
consider the normalisation of the perturbations. Specifi- 
cally, we need the mapping between the amplitude of the 
perturbation and the separation, or orbital frequency, of 
the binary system. The Newtonian tidal potential de- 
ri ved earlier can be u sed to fix this, as follows. 

IPenner et al.l (|2011l ) express the metric component gu 
in the form gu = — e''[l — HQ{r)Yim], with the function 
Hq outside the star taking the form 

Hoir) = apP22{r/M - 1) + aQQ22(r/M - 1) (11) 

where r is the radial distance from the centre of the star, 
P22 and Q22 are the associated Legendre functions, and 
Op, ag are constants. Far from the star the dominant 
term is 



77o(r»i?)«3ap(^— 



(12) 



In this limit the metric can be related to the Newtonian 
potential according to 

1 



Newton 



9tt 



(13) 



This allows us to relate the asymptotic form of Hq to the 
Newtonian tidal potential in equation ([2]), leading to 



ap 



4tT McompM^ 



(14) 



The numerical calculation can then be normalised by 
considering the value of Hq at the stellar surface: 

Hq{R) = ap[P22{,R/M - 1) + a2Q22WM - 1)] (15) 

wh ere the ratio 02 = ap / ap is related to the Love num ber 
k2 (jFlanagan fc Hindered[200l iHinderer et al.ll2010( ) by: 



a2 



(16) 
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The Love number is calculated numerically as a prop- 
erty purely of the stellar model, and the amplitude of 
the associated perturbation is parameterised by ap, ei- 
ther in terms of the binary separation, or the associated 
gravitational-wave frequency: 



ap = Gtt 



C^Mtotal 



comp J.2 

Jgw 



(17) 



Note that, for given stellar masses the tidal deformation 
scales as the square of the gravitational- wave frequency. 

4. RESULTS 

We have generalised the numerical framework of 
IPenner et al.l |2011) to allow for realistic equations of 
state for both the crust and the star's core . We present 
result s for stellar models that combine the lAkmal et al.l 
(I1998D equation of stat e for th e core fluid with the results 
of IDouchin fc Haensell ([200l in the crust. These mod- 
els are state-of-the-art for this problem, but it should be 
noted that we have not accounted for the (likely) pres- 
ence of nuclear pasta in the inner crust. A sizeable pasta 
region could have significant impact on the results, but 
we do not yet have a sufSciently detailed equation of state 
representing this possibility. 
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Figure 1. Left panel: Mass-radius relation for our sequence of 
stellar models, demons trating cons istency with the observational 
cons tra ints ofPcmores t et al.l I|201CI| ) (upper dashed horizontal line) 
and ISteiner et al. (20l3) (grey region). Right panel; The Love 
number k2 (upper curve) as a function of the stellar compactness 
M/R; cf. , the pure fluid results in Fig. 1 of Hindcrer et al. (2010). 
We also show the relative influence of the crust on the tidal de- 
formability, represented by Afc2/fc2""^ (lower curve); this is similar 



to the results of IPenner et al.l l[2011 
IAMq model considered in Fig.|2]is i 
line. 



20111 V The compactness of the 
indicated by a vertical dashed 



The chosen core equation of state is sufficiently stiff to 
satisfy constraints from observations (cf., the left panel 
of Fig [Ij . It allows for neutron star masses at least as 
large as 2Mq, in agreem ent with the observed 1.97M0 
mass of PSR J1614-2230 UDemorest et al.ll2010[) . It also 
satisfies the radius constraint from X-ray burst sources, 
i.e., that a star with ma ss of lAMc^ should have a radius 
in the range 11-12 km (jSteiner et al.112010 '). The elastic 
properties of the crust do not affect the equilibrium con- 
figurations since we assume the star is relaxed at large 
binary separations. 

Given this equilibrium configuration, we have calcu- 
lated both the Love number fc2, and the fractional differ- 
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Figure 2. The gravitational-wave frequency (in Hz) at failure 
for different locations in the crust of a 1.4Mq neutron star. The 
result, which is obtained by com paring the vo n Mises strain to the 
breaking strain of .Horowitz &: Kadaul l|2009l 'l , corresponds to an 
equal-mass binary. The vertical line near 11.15 km indicates the 
location of neutron drip in the star. 

numbers for an elastic s tar an d the equivalent purely fluid 
star (see lPenner et al.l ()201lD for details). The results are 
shown in the right panel of Fig. [TJ They can be com pared 
to, fir st of all, the fluid star resu lts of [H indcrer et alJ 
(|201Clf ) and secondly the results of Penner et al.i (|2011[ ) 
for the magnitude of the crust effects. Based on these 
comparisons, the present results are not surprising. 

We have also evaluated the von Mises stress associ- 
ated with the tidal perturbation. By comparing the 
result to th e anticipated breaking strain Ubrcak ~ 0.1 
(jHorowitz k Kadau 2009), we can infer when different 
parts of the crust fail during binary inspiral. A typical 
result is shown in Fig. [21 providing the gravitational- wave 
frequency at failure throughout the crust for a 1.4M0 
star (in an equal- mass binary). The result is not trivial, 
owing to both the nonlinear combination of perturbed 
quantities that enter the von Mises stress and the asso- 
ciated angular functions. 

The large variation in the crustal strain implies that 
failure will occur at different stages during inspiral. From 
Fig. [5] we see that the outer crust (roughly up to neutron 
drip, corresponding to r « 11.15 km in Fig. [2]) fails fairly 
uniformly when /gw ~ 200 Hz. Meanwhile, failure of 
the bulk of the inner crust requires /gw ~ 600 — 800 Hz, 
a factor of two or so below the ISCO frequency (see equa- 
tion (fTO]) ). However, there are also macroscopic regions 
in the inner crust that will not fail before merger. These 
results improve upon the estimates from section 2, show- 
ing the rich structure of the realistic calculation, with 
failure occurring at different stages at different depths. 

5. IMPLICATIONS 

What happens when the crust fails? Again, key in- 
sights are provide d by the m olecular dynamics simu- 
lations of Horowit z fc Kadaul ([2009) ■ The indications 
are that when the critical strain is reached, there is a 
catastrophic failure, with energy released throughout the 
strained volume, rather than the formation of a lower- 
dimensionality crack. What happens next is less clear. 
Two extreme scenarios can be envisaged: the relieved 
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strain is dissipated locally as heat, or, the strain is 
converted into phonons/seismic waves and transported 
throughout the star prior to dissipation. An interesting 
question is whether or not the relieved strain is capable of 
melting the crust. From the estimates in Section 2 we see 
that the strain energy may be at the level of 10^® erg. 
However, it is easy argue that even the release of this 
level of energy will not have significant impact on the 
crust. To demonstrate this, let us assume that the first 
of the two extreme cases described above applies, with 
all strain energy being dissipated locally. 

Ignoring the temperature dependence of the shear 
modulus (which would be irrelevant for typical neu- 
tron star tempe ratures anyway), the crust melts when 
r < 173, where (|Farouki fc Hamaguchil [19931) 



Here, e is the electron charge, Z is the charge per ion and 

1/3 

a = (S/iTTrii) ' , with rii the ion number density, is the 
average ion spacing. For a typical temperature of 10* K, 
we have F ~ 10* in the inner crust (where Z ^ 50 and 
the number of nucleons per ion. A, is several hundred). 

We want to compare the melting temperature obtained 
from Eq. (|T8)) to the temperature reached after releasing 
the tidal strain. For this we need the shear mod ulus, 
which is approximated by (jStrohmaver et allllQQlD 

/iwO.ln, (19) 

a 

The energy per unit volume corresponding to a strain u 
is then 

-Entrain « 10"'a(^)' (20) 

We can turn this into a temperature via the heat ca- 
pacity, Tstrain « Estrain/cy, making use of the est i- 
mate of the heat capacity of ions from Ivan 
Cv — Cv/ui — aks where a is of order unity (for 10* K 
we have a w 5). Thus, we see that the released heat will 
melt the crust if 

0.5(1)'^' (21) 

This exceeds th e br eaking strain estimates of 
IHorowitz fc Kadaul ()2009) . indicating that the re- 
leased heat does not melt the crust during binary 
inspiral, even when all the strain energy goes into local 
heating. 

A more detailed calculation leads to the results in Fig- 
ure [21 which compares the melting temperature to that 
reached after releasing the strain locally. We also com- 
pare to the critical temperature for (singlet) neutron su- 
perfluidity. The message from this comparison is clear; 
the heat that is released when the crust fails is not suf- 
ficient to break the superfluid pairing, either. This im- 
plies that a realistic numerical treatment of neutron star 
binary inspiral should allow for superfluidity and crust 
elasticity prior to coalescence. 

The crust may not melt, but we are still releasing a 
significant amount of energy when the crust fails. Could 
this have observable consequences? Given the antici- 
pated importance of neutron star binary mergers for fu- 




Figure 3. Various relevant temperatures as functions of baryon 
number density ni, throughout the crust. The upper curve is the 
crustal melting temperature, the lower curve the temperature to 
which the crust would be heated by local dissipation of strain en- 
ergy into heat (cf., the discussion in the text), demonstrating that 
the crust does not melt during inspiral. We also show two mod- 
els for the neutron (singlet) critica l temperature (correspo nding to 
models a (solid) and c (dashed) of lAndersson et al.l l |2003) ). show- 
ing that the local heating is insufficient to break superfluidity. The 
vertical line indicates the location of neutron drip. 

ture gravitational-wave observations, and the likely as- 
sociation of such events with short gamma-ray bursts, a 
precursor signal would undoubtedly be interesting. Ob- 
servational evidencefor such precursors was recently dis- 
cussed by iTroia et al.1 (|2010f ). However, the heating of 
the crust should be irrelevant in this respect. The results 
in Fig. [3| show that the outer crust does not heat signifi- 
cantly. The inner crust heats up, but the associated heat 
wou ld not diffuse to the surface before the binary merges 
(see iBlaes et all ()1989D for discussion) . The possibility 
that the energy is released into seismic waves that gener- 
ate Alfven waves in the magnetosphere, eventually lead- 
ing to a gamma-ra y signal, is more interesting (again, see 
IBlaes et al.l (|1989( )). In principle, such an event could be 
as energetic as the largest observed magn etar flare (the 
27 D ecember 2004 event in SGR 1806-20 (|Palmer et al.l 
l2005f )). assuming that the entire strain energy is trans- 
ferred to the magnetosphere. Such a signal could pos- 
sibly be observable from a distance of 100 Mpc. The 
corresponding gravitational-wave signal would be com- 
fortably detectable by the advanced LIGO/ Virgo detec- 
to r network. H o wever , according to the rate estimates 
of lAbadie et al.l (j2010D there is likely to be only of or- 
der one such event per year of observation. According 
to our numerical results, the most energetic failure event 
would precede the merger itself by a fraction of a sec- 
ond. Observations of such precursors would obviously be 
tremendously interesting, but it seems clear from our re- 
sults that they would have to be associated with rather 
unique events. 

Finally, it remains to consider the impact on the 
gravitational-wave signal. The smallness of the effect 
of the elastic crust on the tidal distortion (right panel 
of Fig. [T]), suggests it would require an unfeasibly high 
signal-to-noise observation for the crustal elasticity to be 
apparent, especially given that finite-size effects them- 
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selves are likely to be only borde rline detectable by th e 
next generation interferometers ()Hinderer et alj 120100 . 
However, as we have shown, even at large separations, 
portions of the crust will fail, with the volume of crust 
pushe d into failur e steadily increasing as the inspiral pro- 
ceeds. 'Kochane3 |l992) suggested that this may lead to 
anomalous frictional damping, which may have a relevant 
secular affect on the phasing. A proper calculation of this 
would take the strain field presented here as its starting 
point, but would incorporate a more realistic treatment 
of crust failure, going beyond the simple breaking strain 
criterion to allow for a regime of plastic deformation and 
dissipation. An accurate study of this effect is clearly 
needed to assess its significance for gravitational-wave 
searches. This would be an important but non-trivial 
extension of the work presented here. 
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